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Ionization by relativistically intense laser pulses of finite duration is considered in the framework 
of strong-field quantum electrodynamics. Our main focus is on the formation of ionization super- 
continua. More specifically, when studying the energy distributions of photoelectrons ionized by 
circularly polarized pulses, we observe the appearance of broad structures lacking the interference 
patterns. These supercontinua extend over hundreds of driving photon energies, thus corresponding 
to high-order nonlinear processes. The corresponding polar-angle distributions show asymmetries 
which are attributed to the radiation pressure experienced by photoelectrons. Moreover, our time 
analysis shows that the electrons comprising the supercontinuum can form pulses of short dura¬ 
tion. While we present the fully numerical results, their interpretation is based on the saddle-point 
approximation for the ionization probability amplitude. 

PACS numbers: 32.80.Rm,32.80.Fb,42.50.Hz 


I. INTRODUCTION 

Interference of probability amplitudes is a fundamen¬ 
tal quantum effect which, among others, manifests itself 
when a strong laser field interacts with matter. In his cel¬ 
ebrated paper on strong-field ionization [lj, Keldysh has 
shown how the probability amplitudes emerging from the 
complex-time saddle points lead to the interference pat¬ 
tern in the angular-energy distribution of photoelectrons. 
This idea, further developed for instance in US (for re¬ 
lated recent investigations, see [1 [b i. has initiated theo¬ 
retical investigations of such quantum processes like the 
above-threshold ionization, the rescattering phenomenon 
or the high-order harmonic generation. Along with im¬ 
portant experimental achievements, it has led to a new 
branch of science called attosecond physics Q, which has 
already found a lot of applications in physics, chemistry, 
biology, and medicine. Some of these achievements are 
described in the review articles (see, e.g., M) or in 
the recent collection of articles devoted to the Keldysh 
theory (l4j . 

The Keldysh approach has been further generalized 
in order to account for the interaction of ionized elec¬ 
trons with their parent ions. This has led to the con¬ 
cept of complex-time trajectories (i.e., trajectories along 
complex time determined by the purely classical New¬ 
ton equations) fl5l - fl7l |. or to the concept of quantum 
complex-time trajectories [l8| which take into account 
the electron wave packet spreading or the quantum dif¬ 
fusion during the complex-time evolution, thus eliminat¬ 
ing problems with the rescattering trajectories present 
in the former approach. Theoretical methods that have 
been developed from the ideas put forward by Keldysh 
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are usually called the Strong-Field Approximation (SFA) 
or the KFR theory (see, also [ijlHoj])- Recently, with the 
development of lasers generating relativistically intense 
and very short pulses, similar interference patterns have 
been investigated in the relativistic strong-field quantum 
electrodynamics [2ll423j . In particular, such diverse phe¬ 
nomena have been studied as the vacuum polarization 
induced Young interference [24j . the Kapitza-Dirac ef¬ 
fect [25|, [26| , interference effects in Compton and Thom¬ 
son scattering (see, e.g., [13, HH), and in laser-modified 
Mott scattering [29|, |3Cjj| , the coherent comb structures 
created bv a finite train of short pulses in the Comp¬ 
ton [31, 32) and Breit-Wheeler [33) processes, which can 
be used for the diagnostics of relativistically intense laser 
pulses HI. 

On the other hand, the lack of interference may result 
in an appearance of supercontinuum (35j which, since 
its demonstration in the early 1970’s, has been the fo¬ 
cus of significant research activities. The supercontin¬ 
uum generation has attracted a lot of attention owing 
to its enormous spectral broadening. Thus, resulting in 
many useful applications, among others, in telecommu¬ 
nication and optical _coherence tomography Util . It has 
been shown in Ref. [37] that the supercontinuum in the 
radiation domain spreading over keV or MeV energy re¬ 
gions can be generated during the Thomson or Compton 
scattering. In light of this result, the question arises: Is 
it possible to generate the supercontinuum in the ioniza¬ 
tion spectrum? In other words: Is it possible to choose 
the parameters of a driving pulse such that the energy 
spectrum of photoelectrons does not change rapidly with 
energy, on the scale of tens or even hundreds laser pho¬ 
ton energies? In general, the answer to this question is 
negative as for even a few-cycle driving pulse the spec¬ 
trum of photoelectrons consists of a sequence of peaks 
separated approximately by the carrier laser frequency 
(which, for monochromatic plane waves, are called the 
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multiphoton peaks). The Keldysh theory shows that such 
a structure arises as the result of interference of at least 
two complex-time probability amplitudes. Thus, in order 
to create the supercontinuum, one should have a domi¬ 
nant complex-time saddle point with an imaginary part 
which, over a broad range of electron final energies, would 
be much smaller than imaginary parts of the remaining 
saddle points. In other words, over a broad domain of 
electron final energies, interference of probability ampli¬ 
tudes should be suppressed. The aim of this paper is to 
show that, indeed, such a situation can happen in the rel¬ 
ativistic ionization of atoms or positive ions by circularly 
polarized, short laser pulses. 

The organization of this paper is as follows. In Sec. HU 
we present the theoretical formulation of relativistic ion¬ 
ization based on the Dirac equation and with the full 
account for the spin degrees of freedom. Specifically, in 
Sec. Ill Al we develop a general approach, which for the 
velocity gauge is elaborated in detail in Sec. IIIBI The 
supplementary Sec. Ill Cl is devoted to the analysis of sad¬ 
dle points arising from our formulation, which further is 
used only for the interpretation of numerical results. In 
Sec. ED we present numerical analysis of ionization of 
He + ions by short, circularly polarized and relativistically 
intense pulses defined in Sec. IIII Al The energy and angu¬ 
lar probability distributions are presented in Sec. IIII Bl 
where also the generation of ionization supercontinua is 
discussed. Sec. IIII Cl is devoted to the further analysis of 
the supercontinuum with the result that it can be shifted 
towards high energies by changing the azimuthal angle 
of emission. In Sec. IIII Dl we study the energy depen¬ 
dence of the phase of probability amplitudes and show in 
Sec. IIII El how the nearly linear dependence of the phase 
on the energy of emitted electrons leads to the time delay 
of the electron wave packets synthesized from a particular 
supercontinuum. In Sec. IIVI we draw some concluding 
remarks. An Appendix n contains supplementary materi¬ 
als concerning the normalization of the laser pulse shape 
functions and physical units. 

Throughout the paper, we keep K = 1. Hence, the fine- 
structure constant equals a = e 2 / (Atteoc) . In numerical 
analysis we use relativistic units (rel. units) such that 
h = m e = c = 1 where m e is the electron rest mass. 
We denote the product of any two four-vectors a^ and 
as a ■ b = a^b^ = a°b° - alb 1 - a 2 b 2 - a 3 b 3 (p = 
0,1,2,3), where the Einstein summation convention is 
used. For the four-vectors we use both the contravariant 
(a 0 , a 1 , a 2 , a 3 ) and the standard (ao,a x ,a y ,a z ) = (ao,a) 
notations. We employ the Feynman notation </l = 7 • 
a = 7 for the contraction with the Dirac matrices 
7 M and use a customary notation u = u' 7 0 , where vJ is 
the Hermitian conjugate of a bispinor u. Finally, we use 
the so-called light-cone variables. Namely, for a given 
space direction determined by a unit vector n (which in 
our paper is the direction of the laser pulse propagation) 
and for an arbitrary four-vector a, we keep the following 
notations: = n ■ a, a~ = a 0 — a\ a + = (a 0 + a^)/ 2 , 

and a 1 = a — a^n. Thus, a - b = a + b~ + a~b + — a 1 - ■ b 3 - 


and d 4 x = da: + dx d 2 ^^. 


II. THEORY AND APPROXIMATIONS 
A. General theory 

The relativistic ionization of one-electron atoms or ions 
is described by the Dirac equation (e = — \e\) 

- e^ v A v [x) - m e c)'F(a:) = 0. (1) 

The electromagnetic potential A(x) is assumed to be of 
the form 

A v (x) = ‘—V(x)5 v o + Ar(x), (2) 

ec 

where V(x) is the binding potential and A^( x) describes 
the laser pulse. The analysis of the time-evolution of the 
system leads to the exact expression for the probability 
amplitude 

Ar = — i J d^e -1 ^ 0 / 0 ^ 4 r f(a;)e^ R (x)4 f i(a:). (3) 

Here, the bispinor \H;(a;) describes the electron bound 
state of energy E 0 and 47 ( 2 ;) is the exact solution of 
Eq. dU) corresponding to the scattering state with the 
incoming spherical waves. 

Although the analytic form of T;(a:) is known for the 
Coulomb potential, the exact solution 47 ( 2 :) can be de¬ 
termined only numerically. Recently, we observe a signifi¬ 
cant progress in solving numerically the relativistic Dirac 
equation (see, e.g., Refs. [38j - |42j ~). Even so, for very large 
intensities of laser pulses available today (of the order of 
10 20 W/cm 2 and larger) and for high energies of photo¬ 
electrons (i.e., a few keV and higher) such solutions are 
not achievable. If, however, the exact scattering state is 
labeled by the asymptotic momentum p for which the 
kinetic energy, y^m^c^J^+Tcp) 2 — m e c 2 , is much larger 
than the ionization potential of the initial bound state, 
m e c 2 — E 0 , then for the final scattering state the Born ex¬ 
pansion with respect to the static potential V ( x ) can be 
applied. This is the essence of the relativistic SFA which, 
in the lowest order, consists in replacing in Eq. © the ex¬ 
act solution 47 ( 2 :) of Eq. © by the corresponding exact 
solution of the simplified Dirac equation 

- eJfi K (x)-m e c)dif\x) = 0 . (4) 

For the relativistic SFA approach to ionization by plane 
wave fields, we refer the reader to [43l - [47| . In this paper, 
however, we focus on relativistic ionization by finite laser 
pulses. 

Let us choose as the solution of (0]), denoted since now 
on by ’Fp^Oc), the one with the well-defined electron mo¬ 
mentum p and spin polarization A = ±. The probabil¬ 
ity amplitude A& ©, which in this case we denote as 



















3 


A(p, A; Ai), becomes 


A{p, A; Ai) = -i J J d 4 xe iq ' x ^°l(x) 

xe4 R (a:)#i(q), (5) 

where 

*i(*) = / (03 eiQ ^i(9)- ( 6 ) 

To shorten the notation, we have introduced q = 
(■ q°,q ) = ( E 0 /c,q ). Note that it is not the four-vector 
as it does not transform properly under the relativistic 
Lorentz transformations. Nevertheless, for the sake of 
space, we shall call it the four-momentum as it is usu¬ 
ally done for the electromagnetic potential A^(x) in the 
Dirac equation m, which is the four-vector only for par¬ 
ticular gauges. In addition, A; labels the spin degrees of 
freedom for the initial state. 

Even after this approximation, the direct numerical 
analysis of Eq. 0 is not possible as, for the laser pulses 
concentrated in a finite focus region, the numerical solu¬ 
tion of Eq. 0 is not available for laser pulse and electron 
parameters mentioned above. For this reason, the so- 
called plane-wave front approximation for the laser beam 
is applied. It accounts for the finite time and space de¬ 
pendence of the laser pulse in the direction of its propa¬ 
gation, however, in the space directions perpendicular to 
the propagation direction it is assumed that the laser held 
extends to infinities. Such an approximation is justified 
if either the laser pulse is not tightly focused or interac¬ 
tion takes place with the highly energetic beams of par¬ 
ticles .atonis or ions in the head-on kinematics (see, e.g., 
Ref. [48|]). Note, that this approximation is commonly 
used in majority of investigations of quantum processes 
in strong laser fields, as only in this case it is possible 
to obtain the exact analytical solution of the Dirac equa¬ 
tion, called the Volkov solution [4f| . 


we have also introduced T p for the pulse duration. The 
two real functions, ), called the shape functions, are 
arbitrary functions with the continuous second deriva¬ 
tives that vanish for <j> < 0 and > 2tt. Introducing the 
relativistically invariant parameter 


Mo| 

fl 5 

m e c 


( 8 ) 


that defines the intensity of the laser field (the relation 
between this parameter and the time-averaged intensity 
of the laser pulse is presented in the Appendix J, we 
rewrite 0 as 


eA(4>) = —m e cp[£\fi{4>) + e 2 / 2 (0)]- (9) 

Hence, the Dirac equation 0 becomes 

{l$ + ITleCjU^/i^) + ^ 2 / 2 (0)] - TOec}^( x) = 0. (10) 

We denote as ( x ) Volkov solution of this equa¬ 
tion for electrons (the superscript (+) means that it is a 
positive-energy solution), where p is the electron asymp¬ 
totic momentum and A = ± labels the spin degrees of 
freedom. Its explicit form can be written as [50 ] 




m e cp 
2p ■ k 


[fi(k ■ x)^t 


+ h{k • M 2 ^]) ( 




( 11 ) 


where 


(x) = p ■ x + 


pk-x 


d 4> 


m e cp 
p ■ k 


( £ i -Ph{4>) 


n-phw)+^f(m+m) 


( 12 ) 


B. Velocity gauge 

In nonrelativistic SFA one usually considers probabil¬ 
ity amplitudes in two gauges: the velocity and length 
gauges. The same is possible in the relativistic case. In 
this paper we consider the velocity gauge, postponing the 
consideration of the length gauge to the near future. 

In order to derive the corresponding formulas let us 
consider the most general form (up to the gauge trans¬ 
formation) of the electromagnetic potential in the plane- 
wave front approximation, 

Ar{x) = A{4>) = A o |ei/i(0) + £ 2 / 2 ( 0 )], (7) 

where <p = k ■ x = k°x~, k = krn = fc°(l,n), fc° = cc/c, 
u) = 2n/T p , £j are two real polarization four-vectors nor¬ 
malized such that £j-£ji = —Sjji and perpendicular to the 
propagation direction of the laser pulse, k ■ £j = 0. Here, 


Here, the photoelectron asymptotic energy equals E p = 
m e c 2 ) 2 whereas the on-mass-shell four-vector 
is p = (p°,p) = (E p /c,p). The Dirac free particle 

bispinors are normalized such that u^u^, = 5w, 
and V is the quantization volume. With this normaliza¬ 
tion the final density of electron states, for a given spin 
degree of freedom, is equal to Vd 3 p/(27r) 3 . 

For our further purposes, we introduced the follow¬ 
ing functions (the explicit form of the ground state wave 
functions v k;(a;) for the hydrogenlike ions and for two spin 
polarizations can be found in the textbook [5l|): 

43,(9) =43^(9), 

4a;A,(9)=4a^ 2*i(9)- (13) 

This allows us to represent the probability amplitude 0 
in the form (note that t jt$ j = I, for j = 1, 2, and 
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tiH\ = °), 


APi A; Ai) 


d 3 q 
(2^)3 


d *xe iS * +)i - x) ~ iq - x M(k-x), 

(14) 


This allows us to write down the probability amplitude 
in the form 


A(p, A; A;) = i rrieCfj, 



V(p, A; Ai), 


(24) 


with 


M(k • x) = i rrieCfj, 

x [/i(fc • x)B ( p \°l (q) + f 2 {k ■ x)Bp^l (g) 

- J^(L/i( fc ' x )] 2 + ^ 2 ( fc ' a? )l 2 ) jB SS i (®)] • ( 15 ) 

Next, we introduce the so-called laser-dressed momentum 

mini 

-P(fi) +£ 2 -P(h))k 

p • k 

+{j %fT -«A 2 > + (16) 

where, for any continuous function F((f>), that vanishes 
outside the interval 0 ^ ^ 27t, we define 

{p) = h r dw) - (i?) 

The dressed momentum fulhlls the equations 

p~ = p~ and p ± =p ± , (18) 



that allow us to write 

S^ + \x) = p + x~ +p~x + — + G p (k°x ~), (19) 

where 


G P {4>) = / del/ 


m e cpL 
p ■ k 


(si -Pif l (</>') - </i)) 


+ S2-p(f2W) - (/ 2 ») + %^(A 2 (^) 


2 p • k 


-(/i 2 > + / 2 2 (</>')-(/ 2 2 )) 


( 20 ) 


Inserting © into © we see that the integrations over 
dcc + and d 2 x^ lead to the conservation conditions 

p~ = q~ and p ± = ( 21 ) 

that permit us to carry out the integration over d 3 g. The 
remaining integration over da;” is performed by applying 
the following Fourier decompositions for j = 1,2 and 
0^0 = k°x~ ^ 27t, 

OO 

[/i(0 )] 3 exp[iG p (^)] = E G^’ 0 ) e” lJV<)i , (22) 

N=—oo 

oo 

[MtfY eMiGM)} = E G^ j) e- iN *. (23) 

N =—oo 


with 


V(p, A; Ai) = E 

iV=^; 

7(1,0) p(l.o) 


2 7ri(p+-q+-)Vfc 0 )/fc 0 _ ^ 


i(p+ — < 7 + — 7Vfc°) 


x \G 


N 


B y ;a(Q) 




m e cp 
2 p ■ n 


[G 


(2,0) 

N 


G 


( 0 , 2 ) 1 jB ( 0 , 0 ) 


N 


pA; Ai 


(Q) (25) 


and 


Q = P+(9°-p°)n- (26) 

Note that in the corresponding nonrelativistic SFA and 
with the dipole approximation applied to the laser field 
the momentum Q in Eq. (l26l) is not shifted by the vec¬ 
tor (q° — p°)n, independently of the gauge choice. This 
term, among others, is responsible for the so-called radi¬ 
ation pressure [54| which has recently been discussed in 
the literature (see, e.g., [42l. [55l - i58l ]). Such a momentum 
shift appears also in the so-called Coulomb-Corrected rel¬ 
ativistic SFA (see, e.g., Eq. (24) in Ref. S3)- 
Eq. (1251) allows us to define the most probable en¬ 
ergy absorbed by photoelectrons during the ionization. 
Namely, the ionization probability distribution is maxi¬ 
mum for such N for which the denominator, p + — q + — 
Nk °, takes the smallest values. This happens for those 
N which are as close as possible to the number 

N e ff = cp+ ~ cq+ = N osc cp+ ~ cq+ . (27) 

LO W L 

Hence, the energy which is most probably transferred 
from the pulse to the atomic or ionic system during the 
ionization equals 


•Etr = N eS 0J = L . (28) 

■L^osc 

The number N e s/N osc estimates how many laser pho¬ 
tons, each carrying the energy wl, have to be absorbed 
for the photoelectron to be detected with the final en¬ 
ergy E p . In numerical illustrations presented in this pa¬ 
per these numbers range from 2000 up to 3000, which 
shows the high degree of nonlinearity of processes being 
considered. Note that both quantities N e s and Et r are 
gauge-invariant [53|. Moreover, N e g is invariant for rela- 
tivistically covariant processes. 

Finally, the total probability of ionization equals 

= J frl^A^)! 2 , (29) 
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whereas its triply-differential distribution takes the form 

A 

(30) 


d 3 P(p, A; Ai) 2 (ro e c) 3 . 2 

= v T7Tvr\P\ ■ \ V (Pi^ X i)\ 


dE p d 2 Qp 


( 27 t) 3 c 


We see that in the last two formulas the quantization vol¬ 
ume V cancels. For the purpose of numerical illustrations 
let us introduce the dimensionless distribution, 


Vkx (p) = a 2 m e c 2 


d 3 P(p, A; Aj) 
d£’pd 2 fl p 


(31) 


which is the probability distribution in the atomic units. 
Later on, we shall label the spin degrees of freedom as 4- 
and t for A or A; equal to — and +, respectively. 

In this Section, we have derived formulas for the 
energy-angular probability distributions that are valid 
only for sufficiently large electron kinetic energies Akin, 
for which the condition: 


-Ekin = V (m e c 2 ) 2 + ( cp) 2 - m e c 2 > m e c 2 - E 0 (32) 


is satisfied. Thus, we shall apply current theory to ioniza¬ 
tion of He + ions, with the ionization potential of roughly 
54eV, and we will analyze the energy-angular probabil¬ 
ity distributions for final electron kinetic energies larger 
than IkeV. It is also commonly assumed that the SFA is 
applicable for sufficiently intense laser fields. Specifically, 
when the ponderomotive energy is larger or comparable 
to the ionization potential. This condition is also very 
well-fulfilled in our numerical analysis. 


C. Saddle-point analysis 


with 


go 

9j 

h 



—m e cp 


-P 
k ■ p ’ 


(m e c^) 2 
2k • p 


3 = 1,2, 


(35) 


Since the integrations over dx + d 2 a;- L lead to the conser¬ 
vation conditions m, we end up with the following ex¬ 
pression for the probability amplitude, 


1 r 2n 

A(p, A;Ai) = ^ J d0e iG (s°’“^)[M(0)] q=g , 

(36) 

where Q is defined in Eq. (12611 . This is the formula which 
suits for the saddle-point analysis. In the following, if it 
does not lead to misunderstandings, we abbreviate the 
function G(go, gi, g 2 , h\ <t>) by G(<p). Let us also note that, 
due to the conservation conditions (1211) , g 0 in (l35l) equals 


go 


p°~q° 
k° ’ 


(37) 


and it depends only on energies of the initial and final 
states. 

Applying now the standard asymptotic procedure for 
the approximate evaluation of integrals, we determine the 
saddle points by solving the equation 


G'(<j>) = 0, 


(38) 


Expressions for the energy-angular probability distri¬ 
bution of photoelectrons, which have been derived in the 
previous Section, are not convenient for the interpreta¬ 
tion of calculated results. On the other hand, a very 
appealing interpretation can be provided by analyzing 
the saddle points of the corresponding integrands, as it 
has been suggested by Keldysh. The latter has been ap¬ 
plied, for instance, in Refs. (4 Ell in investigations of the 
diffraction/interference structures in the probability dis¬ 
tributions. As the formulas presented above do not suit 
for such an analysis, therefore, we have to rewrite the 
expression for the probability amplitude A(p, A; Ai) in a 
different form. For this purpose we present ED in the 
light-cone variables as 

A{p ’ A; Ai) =¥ l # / Wfj dx+d2x± 

x e i(p _ -q~)P 

x (33) 

where M(c/>) is defined by (fl5l) and 

G{go,gi,g2,h-(/))= [ d(/)'[g 0 + 31/1 (</>') +52/2(</>') 

Jo 

+ ( 34 ) 


where ‘ prime’ means the derivative over <j>. This equa¬ 
tion has in general complex solutions. Among them we 
select only those saddle points, denoted by cj> a , for which 
Im G((f> s ) > 0. Hence, we arrive at the approximate ex¬ 
pression for the probability amplitude, 

(39) 

If [W)] 9= is singular at the saddle point, we have 
to apply the so-called singular saddle-point approxima¬ 
tion described, for instance, in Refs. |3l. fl7Ul8l|. However, 
independently of the method applied, the dominant be¬ 
havior of the integral (1361) is determined by the exponent 
e i G Os) which usually, for high energy electrons, decays 
very fast to 0. Moreover, the analysis of such integrals 
shows that in most cases at least two saddle points con¬ 
tribute significantly to the above sum, which results in 
the interference pattern observed in the probability dis¬ 
tributions. At this point, let us emphasize that it is not 
our aim to compare the results predicted by the exact 
formula (l30l) with the ones that follow from the saddle- 
point approximation (13^1) . We treat the expression (l39l) 
only as the appealing interpretative tool for our numeri¬ 
cal analysis. 
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Eq. (l39l) suggests that the interference pattern in ion¬ 
ization is suppressed (i.e., the ionization supercontinuum 
may appear) if there is only one saddle point for which, 
over a broad range of electron energies, ImG(</> s ) is much 
smaller than the corresponding values for the remaining 
saddle points. We shall see below that such a situation 
can indeed take place. Note that conclusions which fol¬ 
low from this interpretation are gauge-independent, as 
only the function M (</>) depends on the chosen gauge. 

III. NUMERICAL ANALYSIS 
A. Pulse shape 

The laser pulse with the sin 2 envelope considered in 
this paper is define as follows. First, we introduce two 
angles, 

Sj = U- i)f, i = 1,2, (40) 

and two functions, 

Fj{$) = F 0 {^,5 j ,x)cos(6 + S j ), (41) 

with 

Fo(4>, Sj,x) = Flo sin 2 sin {N oac ^ + 6j +x) (42) 

for 0 < </> < 2n and 0 otherwise. Next, we define the 
shape functions fj(4>) of the electromagnetic vector po¬ 
tential © as 

W) = ~ [ M'FiW). (43) 

Jo 

Since for N osc > 1 the Fourier decompositions of 
do not contain constant terms, the above definitions of 
) guarantee that also for > 2i r the vector potential 
vanishes. The physical interpretation of the remaining 
parameters is as follows. The angle S determines the 
polarization properties of the laser pulse and, for the cir¬ 
cularly polarized field, we choose <5 = 7 r/ 4 . The carrier 
envelope phase \ is assumed to be 7 t/ 2 . Finally, Nq is the 
normalization-dependent real and positive factor which is 
chosen such that the normalization condition (IA.2I) is ful¬ 
filled. The integer N osc defines the number of cycles in 
the pulse and it is assumed to be equal to 4. The laser 
pulse propagates in the z direction (k = (w/c)n, n = e z ) 
and the real polarization vectors are equal to £i = e x and 
e 2 = e y . 

We consider here the interaction of a laser pulse with 
a positively charged He + ion, which is a one-electron ion 
with Z = 2. Since in our case the pulse is finite, therefore, 
the initial bound state is well-defined. For this reason, 
we do not have to restrict our analysis to one-electron 
ions of very large Z, as it has been done in Refs. 0 - 
113 where the infinite plane-wave field has been consid¬ 
ered. Moreover, as the ion can have relativistic energy, 



-0.5 0 0.5 

-A x (4>)/A s 

x 10~ 5 



£e(0)/£s X 10 5 


Figure 1. (Color online) Trajectories of the tips of the electro¬ 
magnetic vector potential A(0) (upper panel) and the electric 
field vector £((/>) (lower panel) in the relativistic units for the 
laser pulse parameters discussed in detail in Sec. IIII Al All tra¬ 
jectories start from and end up at the origin (0,0). In order 
to show the direction of the time-evolution we mark with the 
filled and open circles the points on these trajectories corre¬ 
sponding to some particular values of </>: red (dark) filled circle 
for rf> = 0.32-7T, cyan (gray) filled circle for (j> = 0.3967T, cyan 
(gray) open circle for (f> = 2n — 0.3967T, and red (dark) open 
circle for 0 = 2zr — 0.327T. The meaning of these particular 
points is discussed in the text. We observe the azimuthal sym¬ 
metry of the electromagnetic potential, p —» n — p mod 27r or 
(. x , y) —> (— x , y), and of the electric field, p —> —<p mod 27t or 
(x,y) (x,-y). 


we can choose the carrier frequency of the laser pulse, wl, 
freely. This is due to the fact that the calculations are 
carried out in the ion reference frame in which, for the 
head-on geometry of the ion and the laser beam, both 
the frequency and the electric-field strength are Doppler- 
upshifted. This aspect has been originally discussed for 
the laser-induced Bethe-Heitler process in Ref. @ (see, 
also the review articles [2l], [13 ). For this reason, we as¬ 
sume that wl = N 0sc w = 20eV and the time-averaged 
intensity is I = 10 20 W/cm 2 (as it has been chosen for 
instance in Ref. 0); in other words, we consider param¬ 
eters which are experimentally available (6^113 ■ 

In Fig. Q] we present trajectories of the tips of vectors 
A(<j>) (upper panel) and £((/)) = —d t A(cj>) (lower panel) 
in the polarization plane xy, which is perpendicular to 
the direction of propagation of the laser pulse n. The 
plots are in the relativistic units, with the units for the 
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E p — m e c 2 (keV) 




6 p /n 



Q p /-k 


Figure 2. (Color online) The ionization probability -P U (P) 
(upper panel) of He + ions for 9 P = 0.487T and ip p = 0 as a 
function of the photoelectron kinetic energy E p — m e c 2 , for 
the laser pulse parameters described in Sec. lIII K\ We observe 
a very broad supercontinuum correlated with the minimum 
of the imaginary part of G(go, gi, g 2 , h\ <j>a) for a particular 
saddle point (lower panel). Im G(go, g\, g2, h; (j> s ) for the re¬ 
maining saddle points are at least two orders of magnitude 
larger. 


Figure 3. (Color online) Upper panel shows the ionization 
probability ( p ) as a function of 9 P for He + ions with en¬ 

ergy E p — m e c 2 = 5565eV and ip p = 0, and for the laser 
pulse parameters described in Sec. IIII Al As the result of the 
radiation pressure, the distribution is shifted towards smaller 
9 P angles compared to the expectations based on the non- 
relativistic theory. This result agrees with the saddle-point 
analysis of the probability amplitude (lower panel), where we 
see that the maximum of the probability amplitude coincides 
with the minimum of lmG(go, gi, g 2 , h; <p s )- 


electromagnetic potential As = m e c/ \e\ and the electric 
field strength £$ = m e c 2 /( |e|Ac) = m^ 3 /l e |i the latter 
is known as the Sauter-Schwinger critical field [63) (Til . 
Hence, in Fig. [T] we present x and y components of 
the vectors —A((j))/As = eA(<j>) / (m e c) and — £(</>)/£§ = 
e£{4>)/(\e\£s). 

The intensity of the laser pulse is considered to be rel¬ 
ativistic if the time-averaged ponderomotive energy (see, 

e-g., 0), 


L r c2 (A-A)-(A).(A) 
2 m e 


(44) 


is comparable to or larger than m e c 2 , or the amplitude of 
the electromagnetic potential is comparable to or larger 
than As- As follows from the upper panel of Fig. [T| for 
the chosen laser pulse parameters the laser pulse intensity 
can be considered as nearly relativistic. 


B. Supercontinuum 

Previously, the main reason to investigate ionization 
of atoms or ions within the relativistic quantum mechan¬ 
ics was to analyze the electron spin effects, not present 
in its nonrelativistic counterpart. It has been shown in 
Ref. 0] that the spin-flipping processes, in which the 
electron initial and final spins both projected on the 


same direction in space (usually chosen as the direction 
of propagation of the laser pulse) are opposite to each 
other, are less probable by roughly two orders of mag¬ 
nitude (i.e., of the order of the fine structure constant 
a) as compared to the ionization processes in which the 
projections of the initial and final electron spins are con¬ 
served. Our numerical analysis confirms these findings 
also for very short pulses. Therefore, in the remaining 
part of this paper we shall analyze the process in which 
the initial and final spins are anti-parallel to the propa¬ 
gation direction of the laser pulse. The second dominant 
process, with the spin projections parallel to the laser 
field propagation direction, only marginally differs from 
the first one. 

Recently, due to the experimental results reported in 
Ref. 0 concerning the effects related to the radiation 
pressure exposed on the atomic or ionic systems by in¬ 
tense laser pulses, we observe an increased interest in the¬ 
oretical investigations of ionization with the relativistic 
effects accounted for (12, (bf] ItI . We shall discuss below 
that the signatures related to the radiation pressure are 
present in the developed above theoretical approach [cf. 
our comments below Eq. (|26D ]. although the main focus 
of our studies is on the formation of a broad ionization 
supercontinuum. 

We consider the probability distribution of electrons, 














P||(p)cos(pp x i cf 4 

Figure 4. (Color online) The ionization probability Vn(p) of 
He + ions for E p — m e c 2 = 5565eV and 9 P = 0.48-7T as the 
function of ip p , for the laser pulse parameters described in 
Sec. HITAl The vertical scale is much smaller than the hori¬ 
zontal one, meaning that the distribution is tightly elongated 
around the x direction. 


'Pli(p), in a given space direction defined by the po¬ 
lar and azimuthal angles, 9 p = 0.487T and ip p = 0, re¬ 
spectively. For photoelectron kinetic energies a little bit 
smaller that IkeV, which still might be in the domain 
of applicability of theoretical methods developed above, 
we observe a typical oscillatory dependence of V^(p) on 
the energy. This can be interpreted as the interference of 
probability amplitudes emerging from at least two sad¬ 
dle points. These oscillations, however, gradually dis¬ 
appear with increasing the photoelectron energy. Such 
that, in the broad part of the spectrum covering few 
hundreds of the laser carrier frequencies wl, we observe 
a smooth behavior of the probability distribution with 
the clearly visible maximum for the electron kinetic en¬ 
ergy E p — m e c 2 ~ 5565eV, as presented in Fig. [2] (upper 
panel). As anticipated above, the interpretation of these 
findings can be based on the analysis of saddle points of 
the probability amplitude (IMl) . Indeed, in the whole do¬ 
main, extending from 3keV up to 9keV, there is only one 
saddle point of imaginary part of G(go, gi, <? 2 , h; <p s ) very 
close to 0, as compared to the remaining saddle points 
for which ImG(go, gi, < 72 , h; <j> a ) are at least two orders of 
magnitude larger. 

Let us analyze further the angular distribution for a 
particular kinetic energy E p — m e c 2 = 5565eV, that ap¬ 
proximately corresponds to the maximum in the super¬ 
continuum presented in Fig. [2] First, we consider the 
polar angle distribution for ip p = 0. The nonrelativis- 
tic theory predicts that for circularly polarized laser field 
ionization occurs predominantly in the plane perpendicu¬ 
lar to the direction of propagation. In our case the max¬ 
imum of the probability distribution should appear for 
6 p = 7t/ 2. However, for sufficiently intense laser pulses, 
due to the so-called radiation pressure, the maximum 
should be shifted towards the direction of propagation, 
i.e., towards smaller angles 6 p . This effect is indeed visi¬ 
ble as presented in Fig. [3] and can be also attributed to 
the minimum of Im G(<?o, gi, 32 , h; <f> s ). 

Although the polar angle distribution of ionization has 
an expected shape if we account for the laser radiation 



0 0.5 1 1.5 2 



Figure 5. (Color online) Shows ImG(jo, gi, < 72 , h; <f> s ) for E p — 
m e c 2 = 5565eV and 9 P = 0.48-7T. We observe that, for a 
given azimuthal angle <p p , there are five saddle points (lower 
panel). Among them only one, let us say <j> 1 , has a very 
small positive imaginary part of G(go, gi, < 72 , h; (f) s ) which is 
observed for ip p = 0.066-7T and tp p = 7r — 0.0667T. These values 
correspond to the real parts of the saddle points Re <j> 1 = 1.687T 
and Re = 0.32-7T, respectively, as shown in the upper panel. 


pressure, the azimuthal angle distribution for relativis¬ 
tic intensities and the high-energy part of the spectrum 
differs from the predictions of the nonrelativistic SFA. 
In the nonrelativistic theory for moderately intense laser 
fields, the ionization occurs with the largest probability 
when the electric field strength of the laser pulse is max¬ 
imum (see, e.g., Ref. 0)- Hence, Fig. Q] would suggest 
that the azimuthal angle distribution should be nearly 
isotropic. Such pattern is not observed in the relativistic 
SFA, as presented in Fig. 2) For this particular elec¬ 
tron final energy the azimuthal distribution is peaked for 
ifip = 0.0667T and ip p = 7r —0.0667T. One can interpret this 
result by performing the saddle-point analysis, presented 
in Fig. [5l In the lower panel we observe that, for each 
of the above azimuthal angles, there is one saddle point 
[we denote it as </>i(y> P )] for which the imaginary part of 
G(go, 51 , 52 , h; 4 >\) is very small. It follows from the upper 
panel of Fig.[5]that these particular azimuthal angles and 
the real parts of the corresponding saddle points can be 
grouped in pairs: (ip p , Re <fii(ip p )) = (0.0667T, 1.687t) and 
( 7 r — 0.0667T, 0.327 t). It is usually assumed that the real 
parts of the saddle points determine the escape time of 
electrons from atoms or ions. The values of the electro¬ 
magnetic vector potential and the electric field strength 
for these two particular phases Re0i(<^ p ) are marked in 
Fig. [T] by red (dark) open and filled circles. The positions 
of these points show that the escape of electrons for rela- 
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Figure 6. (Color online) The same as in Figs. [2] [4] and [3] 
but respectively for: 6 P = 0.48-7T and ip p = 1.87T (top panel), 
E p — m e c 2 = 11750eV and 6 P = 0.48-7T (middle panel), and 
E p — m e c 2 = 11750eV and i p p — 1.87T (bottom panel). 


tivistic intensities takes place not when the electric field 
strength is maximum (as it is in the nonrelativistic SFA), 
but for times when the pulse ramps on and off. This hap¬ 
pens at least for ions with not very large Z. 


C. Azimuthal angle dependence 




<d P / 7r 

Figure 7. (Color online) The same as in Fig. O but for E p — 
m e c 2 = 11750eV and 8 P = 0.487T. 


that now the supercontinuum extends from 8keV up to 
16keV with the maximum for E p — m e c 2 ~ 11750eV. By 
fixing now the photoelectron kinetic energy, E p — m e c 2 = 
11750eV and the polar angle, 9 p = 0.48-7T, in the middle 
panel of Fig. [6] we plot the azimuthal angle distribution 
of ionization probability. As we see, the sidelobes have 
been turned downwards, but still the symmetry x —> —x 
is preserved. As in the previous case, we also see the 
effect of radiation pressure as the maximum of the polar 
angle distribution is shifted to the left (the bottom panel 
of Fig. 0. 

In Fig. 0 we present a similar analysis of the imagi¬ 
nary part of G(go, gi, g2, h', <f> s ) as the function of Rec/> s 
(upper panel) and ip p (lower panel). The minima of this 
function are for (ip p , Re(/>i(<p p )) = (1.27T, 0.3967 t) and 
(1.87T, 1 . 6047 t), that correspond to cyan (gray) filled and 
open circles in Fig. [T] 


The width and the position of maximum of the ioniza¬ 
tion supercontinuum can be controlled by changing the 
azimuthal angle. It appears that, for the considered laser 
pulse shape and the azimuthal angle from the upper part 
of the (x,y) plane, i.e., for 0 < tp p < ir, the supercontin¬ 
uum can be shifted towards lower energies. As a result, 
the interference pattern starts building up. If, however, 
we choose <p p from the interval (7 t,27t), the smooth su¬ 
percontinuum appears for larger energies. These findings 
can be also explain by applying the saddle-point analysis. 

In order to illustrate this effect let us consider the en¬ 
ergy spectrum of photoelectrons ejected in the direction 
given by spherical angles 9 p = 0.487T and (p p = 1.87T. This 
spectrum is presented in the top panel of Fig. [6] We see 


D. Global phase of ionization amplitude 

The energy-angular distribution of photoelectrons dis¬ 
cussed above depends only on the modulus squared of the 
probability amplitude, \A{p, A; A;)| 2 . However, the prob¬ 
ability amplitude, as a complex function, is also charac¬ 
terized by the phase, A; A;), 

A{p, A; A;) = exp[i<F_ 4 (p, A; A ; )] | A{p, A; A ; )|, (45) 

$aO, A; Ai) = arg[A(p, A; A;)], (46) 

defined up to a constant term not affecting any physically 
observable quantities. In particular, this constant term 
disappears in the derivative of the phase over the 
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(Ep - m e c 2 )/co L 

Figure 8. (Color online) The derivative of the phase 4 >.a cal¬ 
culated over the electron energy, Eq. (USD , for two directions 
of emission denoted in each panel. For these angles, the en¬ 
ergy distributions are identical and are presented in the upper 
panel of Fig. [6] 


electron energy. This derivative is presented in Fig. [8] for 
two directions of the electron momentum p, for which the 
azimuthal angular distribution shown in the middle panel 
of Fig. [G]is maximum, i.e., for the angles ( 9 p , <p p ) equal to 
(0.487T, 1.87 t) and (0.487T, 1.27r). We see that these deriva¬ 
tives are nearly constant, over a very broad range of elec¬ 
tron kinetic energies. This means that the phase $_4 can 
be well-represented for this supercontinuum energy re¬ 
gion by the linear dependence, 


®a(p, A; A;) 


4 ) o(u.p, A; Ai) T 4 > i(?ip, A; Ai )E P , (47) 


where n p = p/\p\. 

In order to elucidate what are the consequences of this 
approximate formula, we mention the effects originating 
from the frequency-dependent phase, similar to (S3), in 
classical and quantum optics. It is known that the space- 
and time-dependence of radiation pulses follow from their 
frequency distributions (see, e.g., [371 HH). For radia¬ 
tion pulses, the constant term $o does not influence the 
angular-energy and space-time distributions. On con¬ 
trary, the linear term leads to the time-delay of pulses; 
i.e., the larger the larger time-delay is observed. This 
suggests that, as long as Fig.[5]is concerned, the electron 
wave packet released during ionization in the direction 
{0p,(p p ) = (0.487T, 1.87t) is delayed with respect to the 
one propagating in the direction (9 p , tp p ) = (0.487T, 1.27r). 
This is what the above saddle-point analysis predicts, and 
we are going to address this problem in the next Section. 


The probability amplitude for the electron to be found 
at the point x equals 


4'[x,A;A i |^] = 


l m e c _ 


Vd 3 p 


ip-x, 


u{ P \Ap^i^i)^(p)- 


(48) 

Here, it is assumed that the probability amplitude is cal¬ 
culated for such space-time points x, for which the action 
of the laser pulse is over, i.e., A(k ■ x) = 0. The so-called 
filter function E(p) selects only such final electron mo¬ 
menta, from which we are going to build up the wave 
packet. For our purpose, we select only those electrons 
that move in a given space direction no and have ener¬ 
gies from a chosen range. To be more specific, we assume 
that 


F(p) = 0(U max — E p + m e c 2 ) 

x 0(E P - m e c 2 - £ min )(5 (2) (Qp - (49) 


where 9(E) is the Heaviside step function, E m \ n = 7keV 
and Umax = 18keV [cf. the top panel of Fig. [G] . Hence, 
the electron wave packet propagating in the space direc¬ 
tion x = driQ (where d is the distance from the nucleus), 
and with initial and final spins anti-parallel to the laser 
pulse propagation direction, is proportional to 




(‘■Emax'l'Tfl e C 


' -^min^TTleC^ 


c\F p- iE r> t+ 'Ap \ d „(+) 

a£j P e U |pko,i 


Xy/E^\p\A(\p\n 0 ,l-,l), (50) 


where we have neglected all irrelevant prefactors. Fi¬ 
nally, the time-distance probability distribution of these 
selected electrons is equal to 

'Piit, d) = [A^i, d )] d )• ( 51 ) 

Although the last definition applies to positive times t 
and distances d such that k - x = u>t— (io/cjdn ■ n o > 27r, 
nevertheless, by extrapolating its validity to d = 0 (where 
the nucleus is located) we can roughly estimate the es¬ 
cape time from the ion of the selected electrons, as it 
is presented in Fig. [9] We can judge from the posi¬ 
tions of maxima that the escape times for the selected 
group of photoelectrons are around tot = 1.4647T for the 
direction ( 9 p ,ip p ) = (0.487T, 1.87r) and Lot = 0.3237T for 
the direction (9 p ,ip p ) = (0.487T, 1.27t). This agrees rea¬ 
sonably well with the predictions based on the saddle- 
point analysis, where we have estimated these times to 
be uot — Re = 1.6047T and uot = Re^i = 0.3967T, re¬ 
spectively. 

The width of the distributions presented in Fig. [9] are 
not larger than 1CH 2 . This means that, at the time of 
escape, electrons from these supercontinua are created in 
very short pulses lasting for 

9 1 10 -2 7T 17 

At « HT 2 7 t- = —— » 10 -17 s. 52 

uj 5 eV / ri 









11 




0.3 0.31 0.32 0.33 0.34 0.35 


uit/n 


Figure 9. (Color online) The normalized probability distribu¬ 
tion m for electrons escaping in the directions defined by 
the polar and azimuthal angles as denoted in each panel, and 
for the remaining parameters same as in Figs. [6] and [8] 
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Figure 10. (Color online) The same as in Fig. [9] but for the 
distance d = 5000ao from the nucleus, where a o is the Bohr 
radius. 



Then, during the time-evolution, the electron wave pack¬ 
ets spread. However, even at the distance 5000ao from 
the ion their time durations are of the order of 10 -15 s, as 
it is shown in Fig. [10] As expected, the pulses presented 
there are delayed with respect to each other by the time 
approximately equal to 1.157 t/w. This quite well agrees 
with predictions of the Keldysh theory. 


IV. CONCLUSIONS 


We have developed the theoretical formulation for rela¬ 
tivistic ionization driven by an arbitrary laser pulse. This 
has been done within the lowest order Born approxima¬ 
tion with respect to the atomic potential. As an illus¬ 
tration, we have studied ionization of He + by a short, 
relativistically intense, circularly polarized pulse. As our 
analysis shows, it is possible to adjust parameters of the 
pulse and the target system such that high-order nonlin¬ 
ear processes result in appearance of a broad supercon¬ 
tinuum in the energy spectrum of photoelectrons. 

We have further analyzed the properties of supercon¬ 
tinuum electrons. We have checked that they are pre¬ 
dominantly ionized through the spin-conserved processes. 
We have also observed that, in contrast to the non- 
relativistic ionization by circularly polarized pulses, the 
polar-angle distributions of photoelectrons are asymmet¬ 
ric. This has been ascribed to the radiation pressure ex¬ 
perienced by electrons. As we have argued, in the energy 
region of supercontinuum, the total phase of the proba¬ 
bility amplitude of ionization can be approximated as a 
linear function of the photoelectron energy. Therefore, 
the electron pulses composed out of this part of the spec¬ 
trum can be delayed with respect to each other. More 
importantly, the supercontinuum electrons can form very 
short pulses. Despite the fact that they are spreading 
in time, we have demonstrated that these pulses remain 
fairly short. 

We have shown numerically that relativistic ionization 
can lead to the formation of electron supercontinua. This 
is particularly interesting in the context of designing new 
sources of electron pulses. While our analysis is based on 
purely numerical treatment, we have also performed an 
analysis of probability amplitudes based on approxima¬ 
tion via the saddle-point method. The latter has shown, 
for instance, that relativistic ionization can occur with 
significant probabilities at the pedestal of the driving 
pulse. Not, like in nonrelativistic ionization, at the pulse 
maximum. Other predictions of the saddle-point approx¬ 
imation have agreed well with our fully numerical results. 
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Appendix: Remarks on physical units 


The time-averaged intensity of a pulse described in the 
plane-wave front approximation equals 

r . /"///2\ , ^ (u777l e c) 2 / A i \ 

-'—((/1 ) + \/2 )) £ 0 C g ^ ’ (^- 1 ) 

e A 

where ui is related to the pulse duration T p such that 
T p = 27r/w. If the number of cycles is A osc then the car¬ 
rier frequency is wl = N osc lu , and in our analysis we keep 
it fixed. In order to have the time-averaged intensity in¬ 
dependent of the number of cycles in the pulse we choose 
the normalization of the shape functions such that 


(if > + </2 2 ) = 


and, consequently, 


1 (w L m e c) 2 2 
1 = 2 £ ° C ^ ’ 


We can rewrite this relation as 


1 = 


Ael ^ 


TO e C 2 




(A.2) 


(A.3) 


(A.4) 


where we have restored the Planck constant h , and where 
I le i is the relativistic unit of intensity, 

J re i = « 2.324 x 10 29 W/cm 2 . (A.5) 

87raA^. 87raA 3 7 v ; 

Here, Ac = h/m e c is the Compton wavelength divided 
by 27 t. On the other hand, the nonrelativistic unit of 
intensity is equal to 

7 nrel = a 6 I ie i « 3.51 x 10 16 W/cm 2 . (A.6) 

Thus, 



where a 2 m e c 2 ~ 27.21eV is the nonrelativistic unit of 
energy and 


Mnrel — 


Mo| 

am e c 


(A.8) 


Above, am e c = h/ag is the nonrelativistic unit of mo¬ 
mentum and ao is the Bohr radius. 

Eqs. (IA.4I) and (IA.7I) allow one to relate the parameter 
jx in © to the time-averaged intensity of the laser pulse 
given in the standard units of W/cm". 
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